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Overview 



> Problem Description 

> Computational Model Description 

> Results 

• Self-Pressurization of MHTB Tank 

• Cooling of MHTB Tank using Spray 

> Conclusions 

> Future Work 
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Introduction 



• Affordable and reliable cryogenic storage for use in propellant systems is 
essential to meeting NASA’s future exploration goals. 

• Cryogen mass loss occurs when heat leaks into the tank from the 
surroundings. 

• Heat is carried to the interface by natural convection currents -> evaporation 
-> vapor compression^ rise in tank pressure. 

• Pressure control is necessary to keep tank pressure within design limits 
(venting or active control) 

• Predicting self-pressurization and depressurization rates is important for 
designing future tanks and pressure control systems. 
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Problem Description: MHTB Self-Pressurization and 
Spray Bar TVS Ground-Based Experiment 


Tank Internal volume 37.5 m3 

Cylindrical midsection with: 

height = 3.05 m 

diameter = 3.05 m 

2:1 elliptical end caps 

Tank is enclosed in a vacuum shroud 

4 spray bar tubes attached to center tube heat 
exchanger 

NASA TM-2 12926, 2003 

Goal of this work is to simulate first 
self-pressurization and then cooling 
of the tank via spraying cold liquid 
in to the vapor using ANSYS Fluent 
Lagrangian Spray model combined 
with in-house developed UDFs 


Vacuum Chamber 






MHTB 

Support 

Structure 
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Problem Description 


2D axisymmetric 



26.94 


24.92 


Self-pressurization simulation performed 
on 2D-axisymmetric grid. 

Spray Bar Mixing simulation will use 3D 
90° sector grid. 

Spray-Bar/Heat Exchanger assembly is 
approximated as lying along centerline 


22.90 


20.87 


X 


T[K] 

28.97 


Modeling Approach 



Before starting spray run, 2D-axi 
results interpolated to 3D grid and 
self-pressurization continued for a 
short time to ensure no problems 


3D 90° sector 


Mass Transfer [kg/m A 3-s] 
— 0.25 
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Continuity: 


Computational Model Description: 
Equations Solved 

^ + v(/?v) = 0 
dt v ; 



Momentum: 


— (pp)+ V(pvv) = -Vp + v[p eff (Vv + Vv r )]+ pg + F vol 


Energy: 


| {pE)+V(v(pE + p))= v(k eff VT)+ S h 


2 

Volume of Fluid (VOF) model: £ 

Energy and Temperature are defined as mass average scalars: E = — 


Properties: p = £«/>,, p,„ = K s = £«/#, 

<7=1 r/=l <7=1 

Continuity of Volume Fraction of the q-th phase: 


q=l 


1 


a 


A,L* 


kpJ+V-(« 9 p„vJ=5 




Continuum Surface Force (Brackbill et al.l: F ml = £er 

where /%=V-« ^ 2^'+^ 

kg 

Interfacial mass transfer per unit volume: S a = m. • A . 

“ ’| m -sec 

A. = Vcr , is an interfacial area density in 1/m, lh ; is a mass flux vector in kg/(m 2 -sec). 
where or is a volume fraction of the primary phase 
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Sch rage’s Relation : 


Computational Model Description: 
Equations Solved 



2ff Y M \ 1/2 ( Pi PA 

2 -a) \2nRj \Tp 2 Tp 2 ) 


where cr - accommodation coefficient 



M - molar mass of hydrogen 
R - universal gas constant (8.3 14472 J/molK) 

P t and P v - interfacial and vapor pressures, Pa 

T i and T v - interfacial and vapor temperatures, K (assumed that T L = T v = T sat at the interface) 


Particle Energy Equation: 


UDFs used: 


dT p 

tTlorj C-VJ " 

v v v dt 


h ■ A p (T m Tpj Lirip 


VOF (DEFINE_MASS_TRANSFER) 

♦♦♦ Calculate mass transfer using Schrage relation and supply it to Fluent for phase 
interaction at the interface 

Lagrangian spray (DEFINE_DPM_SCALAR_UPDATE, DEFINE_SOURCE) 

♦♦♦ Perform particle tracking in the vapor, remove particles from the vapor domain when 
they reach the interface and add their contributions to the liquid through source terms. 

♦♦♦ Define sources for the spray bar liquid jets. 

♦♦♦ Model heat and mass transfer between particles (droplets) and vapor 
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Computational Model Description: 
Numerical Methods 



> Simulations performed using ANSYS Fluent version 16 

>2D-Axisymmetric formulation for self-pressurization; 3D 90° sector for Spray cooling study 

> Compressible ideal gas 

> Laminar or k-oo SST turbulence model of Menter et. al 

> Temperature dependent properties for vapor and liquid viscosity and thermal conductivity; 
vapor specific heat 

>lnterfacial mass transfer and mass transfer between droplets and continuous phase (vapor) is 
modeled based on Schrage's relation via user's subroutine 

> Surface tension effects via Continuum Surface Force method of Brackbill et al. 

> Contact angle for hydrogen 0° 


> Second Order Upwind scheme was used for discretization of the Turbulence, Energy and 
Momentum equations (cell values) 

> PISO scheme was used for the Pressure-Velocity coupling (cell values) 

> Least Squares Cell Based scheme was used for the gradient calculations (face values) 

> Body Force Weighted scheme was used for the Pressure interpolation (face values) 

> Point Implicit (Gauss-Seidel) linear equation solver with Algebraic Multi-Grid (AMG) method 
was used for solving linearized systems of equations 

> Bounded Second Order Implicit temporal discretization was used with implicit VOF model; First 
Order Implicit scheme was used with explicit VOF model 
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Computational Model Description: 

Grid 


3D 90° sector (184400 hex cells) 


First row: 0.0094 m 
Growth F actor: 1 . 05 
Number of Rows: 17 


First row: 0.002 m 
Growth Factor: 1.2 
Number of Rows: 14 


First row: 0.002 m 
Growth Factor: 1.2 
Number of Rows: 10 



Number of Intervals 
on the Vapor Side: 136 
Average Interval size: 0.02 m 


Number of Intervals: 151 
Average Interval size: 0.025 m 


Number of Intervals: 64 
Interval size: 0.02 m 




Number of Intervals: 20 
Interval size: 0.12 m 


Number of Intervals 
on the Liquid Side: 97 
Average Interval size: 0.03 m 


On all tank walls y + < 5 
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Computational Model Description: 
Initial Conditions 


Beginning of Self-Pressurization 


Beginning of Spray 



T,K 

33,576 
32.179 
30*982 
23*886 
28.383 
27*092 
25*795 
24.498 
23.202 
21.905 
20.608 




Self-Pressurization 

> T field from experimental data 

> Velocity = 0.0 m/s 

> Turbulent Kinetic Energy (in a turbulent run) = 
1.0e-06 m 2 /s 2 

> Specific Dissipation Rate = 100 1/s 

> Interface initialized at 50% or 90% liquid fill level 

> 2D axisymmetric 



T_vapor, K 
52.01 

44.23 

36.44 

28.66 

20.88 

X 



Spray 

> T, V, vof fields from the end of self-press 
simulation interpolated on to 3D 90 degree 
sector mesh 
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Computational Model Description: 
Boundary Conditions 
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Computational Model Description: 
Injection Setup 



^"Injection type: plain-orifice atomizer with 4 
particle streams per injection 

>lnert particle (coupling with continuous phase 
for mass transfer done in the UDF) 

>Standard parcel release method (releases one 
parcel per injection stream, calculates number 
of particles based on the mass flow rate of the 
particle stream) 

> Injection material is liquid hydrogen with 
constant properties at T=21.088K 

>Two-way coupling with continuous phase; 
unsteady particle tracking with flow time step 

> Particle breakup model 

>Spherical drag law model 

>Variable flow rate based on experimental data 



> Injector inner diameter = 0.001702 m 


>Orifice length =0.000711 m 


>Turbulent dispersion of particles: Discrete 
Random Walk model 
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CFD Results: 

MHTB Tank Self-Pressurization 


Pressure, Pa 
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CFD Results: MHTB Tank Self-Pressurization - 
Accommodation Coefficient Effect 



-5.0E-05 


5000 10000 15000 20000 25000 30000 35000 40000 45000 50000 

time after beginning of self-pressurization, seconds 



Medium Grid: 9,246 cells 
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Implicit YOF 


Explicit YOF 


CFD Results: MHTB Tank Self-Pressurization 
Accommodation Coefficient Effect 
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CFD Results: 

Cooling of MHTB tank using Spray 
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CFD Results: Cooling of MHTB Tank using Spray - 



0 5 10 15 20 25 30 35 40 

time after beginning of spray, seconds 




Mass Transfer Rate at the interface, kg/s 
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CFD Results: Cooling of MHTB Tank using Spray - 
Results of the Turbulent VOF Model 
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CFD Results: Cooling of MHTB Tank using Spray - 
Results of the Turbulent VOF Model 


T_ Vapor (K) 


52 02 


44.24 


3645 


2867 


20 88 





Time [s] = 1.0 Time [s] = 15.0 Time [s] = 34.0 


T [K] 

21.88 


21 63 


21 38 


21.13 


20 88 
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CFD Results: Cooling of MHTB Tank using Spray - 
Results of the Turbulent VOF Model 



T_arop [Kj T_drop [K] T_drop [K] 

26.00 26.88 27.75 28.63 29.50 26.00 26.88 27.75 28.63 29.50 26.00 26.88 27.75 28.63 29.50 



Droplet and temperatures at the center plane of injections 
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CFD Results: Cooling of MHTB Tank using Spray - 
Results of the Turbulent VOF Model 


20.88 

28.66 

T_drop [K] 

36.45 44.24 

52.02 



1 


21.15 

20.88 

28.66 

38.94 

T_vapor [K] 

36.45 44.24 

52.02 

r 


1 1 



20.88 

21.25 

20.88 


28.66 


T_drop [K] 
36.45 


44.24 


28.66 


37.24 
T_vapor [K] 

36.45 44.24 


52.02 


52.02 


20.88 

21.31 

20.88 


28.66 


28.66 


T_drop [K] 
36.45 


33.17 

T_vapor [K] 
36.45 



44.24 


44.24 


52.02 


52.02 



Droplet and vapor temperatures and streamlines at the center 

plane of injections 
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CFD Results: Cooling of MHTB Tank using Spray - 
Results of the Turbulent VOF Model 


T_vapor [K] T_vapor [K] 

26.00 26.88 27.75 28.63 29.50 26.00 26.88 27.75 28.63 29.50 


T_vapor [K] 

26.00 26.88 27.75 28.63 


T_drop [K] T_drop [K] T_drop [K] 

26.00 26.88 27.75 28.63 29.50 26.00 26.88 27.75 28.63 29.50 26.00 26.88 27.75 28.63 


-l r- 



Temperature at the horizontal plane in the vapor 



NCSER 


FLUIDS AND COMBUSTION 


NATIONAL 

CENTER 

FOR 

SPACE 

EXPLORATION 

RESEARCH 


CFD Results: Cooling of MHTB Tank using Spray - 
Results of the Turbulent VOF Model 



MHTB Experiment 


CFD 




Location of the spray bar relative to the temperature 

measurement rake 
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CFD Results: Cooling of MHTB Tank using Spray - 
Results of the Turbulent VOF Model 





48.88 
46.08 
43.28 
40.48 
37.68 

34.88 
32.07 
29.27 
26.47 
23.67 
20.87 



Laminar 
34 seconds 


48.88 
46.08 
43.28 
40.48 
37.68 

34.88 
32.07 
29.27 
26.4 7 
23.67 
20.87 



10 20 30 

time after beginning of spray, seconds 


Turbulent 
34 seconds 


Experiment 
CFD Laminar 
CFD Turbulent 


TD1 


40 
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CFD Results: Cooling of MHTB Tank using Spray - 
Results of the Turbulent VOF Model 



T,K 

1 48.88 
46.08 
43.28 
40.48 
37.68 
34.88 
32.07 
29.27 
26.47 

■ 23.67 
20.87 




Turbulent 
34 seconds 


Laminar 
34 seconds 



48.88 
46.08 
43.28 
40.48 
37.68 

34.88 
32.07 
29.27 
26.47 
23.67 
20.87 


30 
29 
28 
27 
26 
25 
24 
23 
22 

0 10 20 30 40 

time after beginning of spray, seconds 
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CFD Results: Cooling of MHTB Tank using Spray - 
Results of the Turbulent VOF Model 




48.88 
46.08 
43.28 
40.48 
37.68 

34.88 
32.07 
29.27 
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23.67 
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Laminar 
34 seconds 



48.88 
46.08 
43.28 
40.48 
37.68 

34.88 
32.07 
29.27 
26.47 
23.67 
20.87 


Turbulent 
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CFD Results: Cooling of MHTB Tank using Spray - 
Results of the Turbulent VOF Model 




21.88 

21.78 

21.68 

21.58 

21.48 

21.38 

21.28 

21.18 

21.08 

20.98 

20.88 



Laminar 
34 seconds 


T_drop, K 



48.88 
46.08 
43.28 
4 0.48 
37.68 

34.88 
32.07 
29.27 
26.47 
23.67 



Turbulent 
34 seconds 



20.87 
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CFD Results: Cooling of MHTB Tank using Spray 
Results of the Turbulent VOF Model 



Time after beginning of spray, seconds 
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CFD Results: Cooling of MHTB Tank using Spray - 
Results of the Turbulent VOF Model 




Time after beginning of spray, seconds 
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Self-Pressurization 


Conclusions 



• A CFD model was developed for modeling self-pressurization of MHTB tank using 
compressible VOF model with custom model for mass transfer between liquid and 
vapor phases. 

• Varying accommodation coefficient from 0.1 to 0.0001 had very little effect on the 
tank pressure predictions. Explicit VOF model allowed use of larger value of 
accommodation coefficient with a need to reduce time step size when the highest 
value was used 


Spray Cooling 


• A CFD model was developed for simulating spray cooling of MHTB tank using 

compressible VOF with Lagrangian Spray model. 

• The laminar and turbulent VOF models resulted in very similar tank pressures that 
agree well with experimental data. 

• The droplets reduce temperature and promote mixing in the vapor region via heat and 
mass exchange during spray. Temperature of the droplets increases when they travel 
in the vapor towards the interface. Passage of the droplets creates a hot spot in the 
areas of higher droplet concentration in the middle of the vapor region. 

• Droplet accommodation coefficient had significant effect on the tank pressure 
decrease with the higher values resulting in faster pressure drops 
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